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Abstract 

Photonic events with large missing energy have been observed in e + e _ collisions at centre-of- 
mass energies of 130, 136, 161 and 172 GeV using the OPAL detector at LEP. Results are presented 
based on search topologies designed to select events with a single photon and missing transverse 
energy or events with a pair of acoplanar photons. In both search topologies, cross-section measure- 
ments are performed within the kinematic acceptance of the selection. These results are compared 
with the expectations from the Standard Model processes e + c~ — » (single-photon) and 

e + e~ — > 1/7777(7) (acoplanar-photons). No evidence is observed for new physics contributions to 
these final states. Upper limits on cr(e+e~ -> XY)-BR(X -> Y7) and cr(e+e" -> XX)-BR 2 (X -> Y7) 
are derived for the case of stable and invisible Y. These limits apply to single and pair production 
of excited neutrinos (X = v* , Y = v), to neutralino production (X = x®, Y = Xi)i and to super- 
symmetric models in which X = x.i an< i Y — G is a light gravitino. For the latter scenario, the 
results of the acoplanar-photons search are used to provide model-dependent lower limits on the 
mass of the lightest neutralino. 
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1 Introduction 



This paper describes a search for photonic events with large missing energy in e + e _ collisions at centre- 
of-mass energies of 130, 136, 161 and 172 GeV. Two different search topologies are used. Cross-section 
measurements and search results from single-photon and acoplanar-photons topologies at y/s = 130- 
136 GeV Q and y/s = 161 GeV have been previously published. Those results have also been used 
to search for excited neutrinos with photonic decays at y/s = 130-136 GeV Q and y/s = 161 GeV 
In the analyses presented in this paper, both the single and acoplanar-photons search techniques 
are based on those previously published by OPAL, but in each case the kinematic acceptance of the 
analysis has been extended to lower energy and more forward angles, and the efficiency has been 
increased by allowing for the possibility of photon conversions. These results supersede our previous 
results. 

The single-photon and acoplanar-photons search topologies presented here are designed to select 
events with one or more photons and significant missing transverse energy, indicating the presence 
of at least one neutrino-like invisible particle which interacts only weakly with matter. Results on 
photonic events without missing energy are presented in a separate paper. || . 

The single-photon search topology is sensitive to neutral events in which there are one or two 
photons and missing energy, which within the Standard Model are expected from the e + e~ — ► ^7(7) 
process. Measurements of single-photon production have been made in e + e~ collisions at the Z° and 
at lower energies 0-8|. Results from centre-of-mass energies significantly above the Z° mass have 
also been reported V.M. The expected visible energies are sufficiently large at present centre-of-mass 
energies that doubly radiative neutrino pair production can lead to two photons being detected; the 
experimental topology therefore includes such cases. 

The acoplanar-photons search topology is designed to select neutral events with two or more 
photons and significant missing transverse energy, which within the Standard Model are expected 
from the e + e~ — ► 1/1/77(7) process. The selection is designed to retain acceptance for events with a 
number of photons, iV 7 , greater than two if the system formed by the three most energetic photons 
shows evidence for significant missing transverse energy. 

These photonic final-state topologies are sensitive to new physics of the type e + e~ — ► XY and 
e + e~ — * XX where X is neutral and decays radiatively (X — * Y7) and Y is stable and only weakly 
interacting. For the general case of massive X and Y this includes conventional supersymmetric pro- 
cesses [[RJ (X = Y = Xi)- I n this context it has been emphasised ]ll| that the radiative branching 
ratio of the x.2 m &y be large. There is particularly good sensitivity for the special case of My ~ that 
applies both to the production of excited neutrinos (X = u*, Y = u) and to supersymmetric models in 
which the lightest supersymmetric particle (LSP) is a light gravitinc[], and Xi 1S the next-to-lightest 
supersymmetric particle (NLSP) which decays to a gravitino and a photon, (X = Xi, Y = G). In this 
case, the branching ratio of this decay of the Xi is naturally large. Such a signature has been discussed 



in [12] and more recently in [13-15] for a no-scale supergravity model and in [16 for a model with 
gauge-mediated supersymmetry breaking; in each case, X1X1 production cross-sections of order 1 pb 
are predicted at these centre-of-mass energies, for Mso ~ 50 GeV. Other types of new physics to which 
these search topologies are sensitive include the production of invisible particles made visible through 
initial-state radiation and the production of an invisible particle in association with a photon. The 
acoplanar-photons search topology also has sensitivity to the production of two particles, one invisible, 
or with an invisible decay mode, and the other decaying into two photons. Such events might arise 
from the production of a Higgs-like particle, S° : e + e~ — > Z°S°, followed by S ^ 77, Z° — ► vv. Results 



1 The mass scale is typically O(keV). 
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for this model searching for the hadronic and leptonic Z° decays appear elsewhere [17 



This paper will first briefly describe the detector, the data sample and the Monte Carlo samples 
used, including a discussion of event generators for e + e~ vV + wy. The event selection for each 
search topology will then be described, followed by cross-section measurements for e + e~ — > 1^7(7) and 
e + e~ — > 1^1777(7) and comparisons with Standard Model expectations. Implications of these results 
on the possibility of new physics processes of the type e + e~ — ► XY or XX, X — ► Y7 will be discussed. 



2 Detector, Data Sample and Monte Carlo Samples 



The OPAL detector is described in detail elsewhere [|18|. The measurements presented here are mainly 
based on the observation of clusters of energy deposited in the lead-glass electromagnetic calorimeters 
(ECAL) consisting of a 9,440 lead glass block array in the barrel (|cos#| < 0.82) with a quasi- 
pointing geometry, and two dome-shaped endcap arrays, each of 1,132 lead-glass blocks with axes 
coaxial with the beam axis covering the polar angle range (0.81 < |cos^| < 0.984). In the overlap 
region, 0.785 < |cos#| < 0.815, and at very forward angles, |cos#| > 0.94 the energy-resolution 
of the ECAL is slightly degraded relative to the nominal resolution. In some cases (where stated) 
these regions have been excluded from the analysis. These calorimeters, together with the gamma- 
catcher calorimeter, the forward calorimeters and the silicon-tungsten calorimeter (SiW), provide a 
fully hermetic electromagnetic calorimeter down to polar angles of 33 mrad. The SiW calorimeter 
covers polar angles down to 24 mrad, however the region around 30 mrad lacks useful calorimetric 
coverage due to the installation, in 1996, of a thick tungsten shield designed to protect the tracking 
chambers from accelerator backgrounds while running at centre-of-mass energies well above the Z° 
resonance. The tracking system, consisting of a silicon microvertex detector, a vertex drift chamber 
(CV) and a large volume jet drift chamber (CJ), is used to reject events with prompt charged particles. 
The silicon microvertex detector consists of two concentric cylindrical layers of silicon microstrip 
arrays, each layer providing both an azimuthal and longitudinal (along the beam direction) coordinate 
measurement. The two layer acceptance covers | cos#| < 0.90 for the 161 and 172 GeV data-set while 
for the data acquired in 1995 at 130 and 136 GeV, the acceptance is limited to | cos 6*| < 0.75. 
Beam-related backgrounds and backgrounds arising from cosmic-ray interactions are rejected using 
information from the time-of-flight system, (TOF), the hadron calorimeter and muon detectors. 

The data used in this analysis were recorded at e + e~ centre-of-mass energies of 130.3, 136.2, 161.3, 
and 172.1 GeV, with integrated luminosities of 2.30, 2.59, 9.89, and 10.28 pb _1 , respectively. The 
integrated luminosities are determined to better than 1% from small-angle Bhabha scattering events 



in the SiW luminosity calorimeter. Triggers [19] based on electromagnetic energy deposits in either 
the barrel or endcap electromagnetic calorimeters, and also on a coincidence of energy in the barrel 
electromagnetic calorimeter and a hit in the TOF system, lead to full trigger efficiency for photonic 
events passing the event selection criteria described below. 

For the expected Standard Model signal process, e + e~ — > vv + wy, the Monte Carlo generators 
NNGG03 §0|, NUNUGPV @ and KORALZ g| were used. Modelling of these backgrounds is 



discussed in more detail in section |2.1| . For other expected Standard Model background processes, 
a number of different generators were used: RADCOR (2j| for e + e~ — * 77(7); BHWIDE |24j and 
TEEGG f2§ for e+e" -» e+e"^); and KORALZ for e+e" -> M + M - (7) and e+e" -» T+T~{i). 



To simulate possible new physics processes of the type e + e — * XY and e + e — * XX where X 



decays to Y7 and Y escapes detection, the SUSYGEN [26 Monte Carlo generator was used to produce 
neutralino pair events of the type e + e~ — > X2X1 an d e + e~ — > X2X2) X2 ~~ * X?7j with isotropic angular 
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distributions for the production and decay of x% and including initial-state radiation. SUSYGEN 
Monte Carlo events were generated at each centre-of-mass energy, for 16-24 points in the kinematically 
accessible region of the (M%, My) plane for which Mx-My > 5 GeV, depending on the centre-of-mass 
energy. For the case My ~ 0, the efficiencies for XY and XX production obtained from these Monte 
Carlo samples are consistent within statistical errors with the efficiencies obtained from the OPAL 
v*v and v*v* excited neutrino Monte Carlo samples respectively |J. All the Monte Carlo samples 
described above were processed through the OPAL detector simulation [27]. 



2.1 Event Generators and Analytical Calculations of e+e — > w + wy 

The present status of event generators for the Standard Model process e + e _ — > vv + 717, n > 1, is 
very unsatisfactory for the centre-of-mass energy region, 130 GeV < y/s < 172 GeV, relevant to the 
analyses presented here. 

For y/s ~ Mz, two event generators, NNGG03 and KORALZ were used extensively for studies of 
e + e~ — > 1/7/7(7) with a demonstrated agreement between the cross-section predictions of better 
than 1%. NNGG03 is designed for e + e~ — ► 1/777(7) at y/s ~ Mz with inclusive exponentiation of soft 
photons and the hard photon matrix element for e + e~ — ► 1^77(7) for the Z exchange diagrams only. 
At higher centre-of-mass energies, it has not been maintained officially by the authors, nor claimed 
to be reliable. The absence of a complete lowest order calculation for e + e~ — > z^z/77 and higher order 
corrections (e + e~ — > 1/1/77(7)) make it necessarily incomplete for e + e~ — ► 1^77(7). The KORALZ 
event generator, primarily designed for e + e~ — > /^ + /i~(7) and e + e~ — > t + t~(^) at y/s « Mz, can 



also generate e + e~ — * i/V + nj, n > 0, using the Yennie-Frautschi-Suura approach [29| to explicitly 
generate an arbitrary number of additional initial state photons. This generator is maintained by 
the authors for y/s 3> Mz, but no specific publications exist yet attesting to its accuracy for either 
1/177(7) or 7/7777(7) in this regime. The NUNUGPV analytical calculation is designed for y/s ~ Mz 
and y/s 3> Mz using the pt dependent structure function approach to estimate e + e~ — > 2/177(7) with 
a claimed accuracy of 1-2% for 150 < y/s < 175 GeV. An event generator based on this calculation 
is also available which includes the emission of an additional photon from each beam. This feature 
is designed to permit estimation of the effect of 1/1777(7) events on the 1/777(7) acceptance. It is not 
intended as an accurate estimate of the 1/7777(7) cross-section. In a previous publication [g], we used 
this feature inappropriately to estimate the expected contribution from e + e~ — > ^77(7). 

Recently, independent calculations have been made of the e + e~ — > 1/F77 (lowest order) cross- 
section using CompHep (3(| (by Ambrosanio |y]), and using HELAS |3^] (by Mrenna f33|). Calcula- 
tions of the e + e~ — > 2/7777(7) cross-section by Mrenna, and by Bain and Pain |34| using GRACE [pq| 
and CompHep have also been made. These are approximately a factor of two lower than the 1/7777(7) 
cross-section predictions we obtain using the NNGG03 and NUNUGPV event generators. The esti- 
mated 1/7777(7) cross-section from KORALZ agrees reasonably well with the independent calculations. 
For e + e~ — ► 1/7/7(7), we have found that the estimated cross-section from KORALZ is lower by about 
10% compared with NUNUGPV within the kinematic acceptance of the single-photon selection, de- 
scribed in section |l[ 

KORALZ is used to estimate the detection efficiency of e + e~ — * 1/777(7) and e + e~ — * 1/7777(7) 
given its more complete treatment of events with multiple photons. Generator studies indicate that 
it also provides a reasonable estimate of the fraction of two photon events (e + e~ — ► 1^77(7)). The 
estimated efficiencies obtained using KORALZ are compared with those obtained using NNGG03 
and NUNUGPV and only small differences are found, indicating that the experimentally measured 
cross-sections are relatively insensitive to the choice of generator. 
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For coherence in the comparisons of data with Monte Carlo, we use KORALZ. In calculating 
upper limits on new processes for the single-photon topology, despite the claimed 1-2% accuracy of 
NUNUGPV, we use the background estimate from KORALZ which is the lower of the two and is 
therefore expected to be conservative. Given the current status of calculations of e + e~ — > 1^77(7), 
where factor of two differences between some cross-section estimates are not yet fully understood, in 
calculating limits on new processes for the acoplanar-photons topology we do not take into account 
the ^7/77(7) background estimate. 



3 Photonic Event Selection 

This section describes the criteria for selecting single-photon and acoplanar-photons events. The kine- 
matic acceptance of each topology is defined in terms of the photon energy, and the photon polar 
angle, 6, defined with respect to the electron beam direction. In order to simultaneously maintain sen- 
sitivity to low photon energies and to retain acceptance at high polar angles, the kinematic acceptance 
of each selection is composed of two (possibly overlapping) parts: 

Single-Photon - One or two photons accompanied by invisible particle(s): 

• At least one photon with xt > 0.05 and | cos 9 \ < 0.82, or, 

• At least one photon with xt > 0.1 and (0.82 < | cosQ \ < 0.966). 

Acoplanar-Photons - Two or more photons accompanied by invisible particle(s): 

• At least two photons with x 1 > 0.05 and 15° < 9 < 165°, or, 

• One photon with > 1.75 GeV and | cos 9 \ < 0.8 and a second photon 
satisfying £ 7 > 1.75 GeV and 15° < 9 < 165° (| cos<9 | < 0.966). 

where the scaled energy, x 7 , is defined as Ey/E\j eajia , and the scaled transverse momentum, xt, is 
defined as E 1 sm0/£'b eam . 

In each of the topologies, it is desirable to retain acceptance for events with an additional photon, 
if the resulting multi-photon system is still consistent with the presence in the event of significant 
missing energy. This reduces the sensitivity of each measurement to the modelling of higher-order 
contributions. 

3.1 Single-Photon Event Selection Description 

Events pass the single-photon selection if they satisfy the criteria listed below. These selection criteria 
are similar to previous OPAL analyses of photonic events with missing energy but have increased 
acceptance and efficiency: 

• Angular acceptance and minimum transverse energy. An event is considered to contain a 
photon candidate if the primary electromagnetic cluster (that with the highest deposited energy 
in the barrel or endcap calorimeters) is in the region 15° < 9 < 165° (|cos#| < 0.966) and 
has a scaled transverse energy, xt, that exceeds 0.1. Events with a primary cluster having 
0.05 < xt < 0.1 and in the barrel region | cos 9 \ < 0.82 are also accepted if they have associated 
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TOF information with good timing, as described in the fourth selection criterion below. Events 
are considered to have more than one photon if additional electromagnetic clusters are found in 
the barrel or endcap calorimeter (| cos# | < 0.984) having deposited energy exceeding 300 MeV. 

• Charged track veto or photon conversion consistency requirements. Events are vetoed 
if there is a charged track with ten or more hits in the central detector, unless the track is 
consistent with arising from a photon conversion. Events having no charged tracks are called 
non-conversion candidates. Photon conversion consistency requires that the primary photon 
candidate cluster be associated with the highest pr track in the event within 100 mrad in both 
azimuthal and polar angle. This track should not be prompt, i.e. the cluster is not accepted as a 
possible photon conversion if there are two or more (out of a maximum possible of four) associated 
silicon microvertex detector hits for photons within its two-layer acceptance, or a minimum 
number of hits in the CV inner axial wires (for photon polar angles beyond the microvertex 
two-layer acceptance) that are associated in azimuth to the above cluster. Events for which 
the primary photon is consistent with a photon conversion are called conversion candidates. 
Events with conversion candidates are rejected if they have at least 2 tracks, reconstructed from 
axial-wire hits in CV, with an opening angle in the transverse plane exceeding 45 degrees. This 
cut removes most of the events having charged tracks which do not arise from a single photon 
conversion. 

• Cluster extent. Only clusters containing more than one ECAL block are considered as photon 
candidates. The primary electromagnetic cluster, combined with any clusters contiguous with 
it, must be consistent with the cluster size and energy sharing of blocks for a photon coming 
from near the interaction point. The cluster size varies in both azimuthal and polar angle extent 
as a function of | cos 9 \ . The cluster extent cuts are parametrized in | cos 6 | accordingly. 

• Forward energy vetoes. The total energy deposited in each of the left and right forward 
calorimeters and in each of the left and right SiW calorimeters must be less than 5 GeV. The 
most energetic gamma-catcher cluster must have an energy of less than 5 GeV. These vetoes 
serve to ensure that photon candidate events are not accompanied by any event activity in the 
forward regions. 

• Muon veto. Events are rejected if there are any muon track segments reconstructed in the barrel 
or endcap muon chambers, or in the barrel, endcap or pole-tip hadron calorimeters. Events are 
also rejected if three or more of the outer eight layers of the barrel hadron calorimeter have 
strips hit in any 45° azimuthal octant. The muon veto is used primarily to remove cosmic ray 
background. 

• Timing measurement in TOF system for low xt and conversion candidates. An 

electromagnetic cluster is said to have an associated TOF hit if it is matched within 50 mrad 
in azimuthal angle by a good quality TOF counter signal produced by the photon converting 
before or in the coil which is located in front of the TOF. A cluster with an associated TOF 
hit has good timing if the measured arrival time of the photon at the TOF is within 20 ns of 
the expected time for a photon originating from the interaction point. For all events with a 
photon conversion candidate in the barrel region | cos 6 \ < 0.82 and for events with a low xt 
(0.05 < xt < 0.1) non-conversion candidate in the barrel region, we require an associated TOF 
hit with good timing. For all other events with the primary photon in the barrel region, if there 
is an associated TOF hit, it must have good timing. 

• Special background vetoes for events with no TOF hit. If the primary non-conversion 
candidate photon has no associated TOF hit, three different background vetoes are applied. The 
first rejects events in which any of the three muon triggers (T^j (barrel and two endcaps) were 
present. This veto removes much of the cosmic ray background. The second looks for a series 
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of calorimeter clusters at the same r and <fi as the primary cluster, but at different z. This veto 
rejects beam halo type backgrounds. The third looks for a series of hits in the outer layers of 
the hadron calorimeter. This veto rejects both cosmic rays and beam related backgrounds. 

Events with a second photon candidate are rejected as background from e + e~ — > 77(7) whilst 
retaining acceptance for events with two photons and missing energy if any of the following criteria 
are satisfied: 

• The total energy of the two clusters exceeds 0.9-y/i. 

• The acoplanarity angle^j of the two clusters is less than 2.5°. 

• The missing momentum vector calculated from the two clusters satisfies | cos# m i ss | > 0.9. 

• A third electromagnetic cluster is detected with deposited energy exceeding 300 MeV. 

• The transverse momentum of the two photon system does not exceed 0.05 of the beam energy. 

• For events with at least one of the two clusters in the region | cos 9 \ > 0.95, the variable 6t is less 
than 0.1, where 6t = (sin #1 + sin #2) I cos [(^i — 4>2)/2] |. This amounts to a stronger acoplanarity 
cut for events with at least one forward photon. 

For the conversion selection, Figure |l| a) shows the maximum of the angular separation, in 6 and 
<j), of the primary photon candidate and the highest track in the event. Overlaid as a histogram 
is the expected distribution from 1/1/7(7) Monte Carlo, normalized to the integrated luminosity of the 
OPAL data. The cut at 100 mrad rejects contributions from cosmic rays. For the non-conversion 
selection, Figure [l] b) shows the difference between the measured TOF timing and that expected for 
a photon from the interaction point, for events passing all selection criteria or failing only the TOF 
timing requirement. The seven events outside the accepted region of ±20ns are rejected as cosmics. 

3.2 Acoplanar-Photons Event Selection Description 

The acoplanar-photons selection has two overlapping regions of kinematic acceptance, in order to 
retain both sensitivity to low-energy photons and acceptance at large | cos 9 \ . These selections are 
based on analyses previously published by OPAL using data collected at centre-of-mass energies of 
130-136 GeV Q. The analysis presented in this paper has increased acceptance and efficiency relative 
to the previous OPAL analyses. The selection criteria are summarized below: 

• Angular acceptance and minimum energy. Events are accepted as candidates if there 
are at least two electromagnetic clusters with scaled energy, x 7 , exceeding 0.05 in in the polar- 
angle region 15° < 6 < 165° (| cos# | < 0.966). In order to retain sensitivity to physics processes 
producing low-energy photons, the minimum energy requirement is relaxed to 1.5 GeV deposited 
energy for photon candidates in the polar-angle region | cos 8 \ < 0.8. This energy deposition 
corresponds to a photon with an energy of about 1.75 GeV ||. These two selections are referred 
to below as the "standard" and "low-energy" selections, respectively. Background vetoes are 
applied differently for the two parts of the selection, as described below. 

2 Defined as 180° minus the opening angle in the transverse plane. 
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Charged track veto or photon conversion consistency requirements. We use selection 
criteria designed to reject events having tracking information consistent with the presence of at 
least one charged particle originating from the interaction point. These criteria are designed 
to retain acceptance for events in which one or both of the photons convert. For the standard 
selection we use hit information from the central jet-chamber, the vertex drift chamber, and the 
silicon microvertex detector (for | cos 9 | < 0.75 (0.9) for data taken at y/s = 130-136 (161, 172) 
GeV). These three detectors form independent estimators for the existence of charged particle 
activity. Events in which charged particle activity is associated in azimuth with both photon 
candidates are rejected unless the signal is from the jet chamber only or from the microvertex 
detector only. Events in which only one photon candidate has azimuthally associated charged 
particle activity are rejected if all (two or three) layers of charged particle detection registered 
activity. To address possible backgrounds from e + e~ —* £ + l~jj, an additional veto requires 
that there be no reconstructed charged track with transverse momentum exceeding 1 GeV, with 
associated hits in the axial layers of the vertex chamber, and separated from each of the photon 
candidates by more than 15°. 

The low energy part of the selection requires that there be no reconstructed charged track in 
the event with 20 or more reconstructed hits in the central jet-chamber. 

Cluster extent. Any photon candidate within the polar angle region | cos 9 \ < 0.75 is required 
to have an angular cluster extent that is less than 250 mrad in both 9 and (f>. 

Forward energy vetoes. The forward vetoes described for the single-photon analysis are 
applied with the same thresholds. 

Muon veto. To suppress backgrounds arising from cosmic-ray muon interactions or beam 
halo muons which can deposit significant energy in the calorimeter, the events must pass the 
muon veto described for the single-photon analysis. Additionally the special background vetoes 
described for the single-photon selection are applied to events in which no TOF information is 
present. 

Timing measurement in TOF system. Requirements on time-of-flight (TOF) information 
are defined separately for the two parts of the kinematic selections. For the low-energy part of 
the selection, we require that the photon in the barrel region has an associated TOF hit with 



good timing (as defined for the single-photon analysis in section ^2). For the standard selection, 
for events in which both of the photon candidates lie within | cos 9 \ < 0.82, at least one of them 
must have an associated TOF hit with good timing. For all events we reject the event if either of 
the photon candidates has an associated TOF hit with bad timing. Finally, if there is a charged 
track associated with a cluster within the polar region | cos 9 \ < 0.82, the requirement of an 
associated TOF hit with good timing is applied. 



Acoplanar photons events can be faked by cosmic-ray and beam-halo events in which a muon grazes the 
electromagnetic calorimeter. Such events can produce large clusters which are split by the clustering 
algorithm to produce two or more clusters. Since it is difficult to model such backgrounds it is useful to 
have a great deal of redundancy in the procedures used to reject these contributions. This redundancy 
is provided by the selection criteria outlined above. Figure [l] c) shows the maximum cluster extent 
for events in which both photons are in the polar-angle region |cos#| < 0.75, where cluster-extent 
cuts are applied. The shaded histogram shows the cluster-extent distribution for real photons coming 
from the interaction point. These come dominantly from e + e~ — > 77 events selected by removing 
the anti-77(7) cuts (described below). The shaded histogram shows the same distribution for events 
failing the TOF cuts and (possibly) the special background vetoes. The cut at 250 mrad is indicated 
by the arrow. Additional suppression of such events, as well as of beam-wall/beam-gas events and 
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instrumental backgrounds in the overlap and endcap regions, is obtained by imposing the following 
requirements: 



• Events are vetoed if there is a reconstructed charged track with at least 20 hits in the jet chamber 
and a \zq\ larger than 50 cm, where zq is the z coordinate of the point of closest approach of the 
track to the beamline in the transverse plane. 

• Events are vetoed if the total number of ECAL clusters having more than 1 GeV of deposited 
energy is larger than five. 

• If both photons have | cos 6 | > 0.75, the opening angle ip of the two-photon system is required 
to satisfy cos^ < 0.98. Otherwise the requirement is that the azimuthal separation of the two 
candidate clusters be greater than 2.5°. 



Finally, background from e + e — > 77(7) is rejected, whilst retaining acceptance for the signal 
topology, if any of the following criteria are satisfied: 



• The total visible energy of the event exceeds 0.95-^/5. 

• The acoplanarity angle of the two highest-energy clusters is less than 2.5°. 

• The missing momentum vector calculated from the two highest-energy photon candidates satis- 
fies I cos # m i ss I > 0.95. 

• The transverse momentum of the two-photon system is less than 0.05Sbeami events having a 
third photon candidate (with EL > 300 MeV) are rejected unless the three photon system is 
significantly aplanar (sum of the three opening angles < 350°) and the transverse momentum of 
the three-photon system exceeds 0. I-beam- 



Figure [T| d) shows the distribution of the acoplanarity angle for events passing all cuts or failing 
one or both of the total-energy cut and the cut on the scaled transverse momentum of the two-photon 
system. The OPAL data is shown as solid points with error bars. Overlaid as a histogram is the 
expected distribution, from e + e~ — > 77(7) Monte Carlo, normalized to the luminosity of the OPAL 
data. The cut at 2.5° is indicated. 



4 Results 



The results of the single-photon and acoplanar-photons selections are given below in sections 4.1 and 
section |4.2| respectively. In each section, the measured cross-sections for the search topology are given 
and compared with Standard Model expectations, and the results of the XY and XX searches are 
then described. Upper limits on a{e + e~ -> XY) • BR(X Y7) and cr(e+e- XX) • BR 2 (X -> Y7) 
are given, respectively. This is done both for the general case of massive X and Y, applicable to 
conventional super symmetric models in which X = X2 an d Y = Xi> and also separately for the special 
case of My ~ 0, which applies both to single and pair production of neutralinos in supersymmetric 
models in which the LSP is a light gravitino and to single and pair production of excited neutrinos. 
These results are used to set limits on the production of excited neutrinos (u*) in a separate paper [pq| . 
For all such limits, the efficiencies were evaluated with the decay length of X set to zero. The efficiencies 
are unaffected if the decay length is much less than the distance from the interaction point to the 
electromagnetic calorimeters (~2 m). 
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For the single-photon and acoplanar-photons analyses, XY and XX Monte Carlo events were 
generated at each centre-of-mass energy for a variety of mass points in the kinematically accessible 
region of the (Mx,My) plane. To set limits for arbitrary Mx and My, the efficiency over the entire 
(Mx, Afy) plane is parametrized using the efficiencies calculated at the generated mass points. For the 
single-photon search topology, the region with Mx + My < Mz is kinematically accessible at yfs ~ Mz 
and strong limits have already been reported |57]. For the acoplanar-photons search topology, we 
restrict the search to Mx values larger than about Mz/2. At lower masses, limits have been reported 
at yfs ~ Mz [38 1 and possible radiative return to the Z followed by Z — > XX would yield very different 
event kinematics than those of the signal Monte Carlo events used for this study. For both search 
topologies, at values of Mx — My < 5 GeV, the estimated efficiency decreases significantly due to 
event kinematics that yield low photon energies. For that reason no limits are set in this region. 



4.1 Single-Photon Results 7(7) + tf T 



After applying the selection criteria of the single-photon selection to the \fs = 130-172 GeV data 
samples, a total of 138 events are selected. The expected non-physics background is 2.3 ± 1.1 events, 
consisting solely of cosmic ray and beam related backgrounds. This non-physics background has 
been estimated from events detected out of time and using a visual scan with looser cuts. The 
expected physics backgrounds from plausible sources, e + e~ — > 77(7), Bhabha events with initial or 
final-state radiation and e + e~ — > fi + fi^ r y and e + e~ — > t + t~j, have been evaluated to be less than 
0.4 events at 95% confidence level (CL) and are therefore considered to be negligible for the cross- 
section measurement. For each of the four centre-of-mass energies, Table 1 shows the number of 
events observed, the number of events expected from the Standard Model process e + e~ — > 2/277(7) 
evaluated using the KORALZ generator, the NNGG03 generator = 130 and 136 GeV) and the 
NUNUGPV generator (y/s = 161 and 172 GeV), and the number of background events expected 
from non-physics processes. The estimated efficiency for selecting e + e~ — > 2/777(7) events within the 
kinematic acceptance of the single-photon selection is also given, as is the corresponding measured 
e + e~ — ► 1/777(7) cross-section within this kinematic acceptance, corrected for detector and selection 
efficiencies, and subtracting the estimated non-physics background. For both the single-photon and 
acoplanar-photons selections, efficiency losses due to detector occupancy range from about (3-5)% at 
the different centre-of-mass energies. Here and elsewhere in this paper, unless otherwise stated, quoted 
efficiencies include these losses and those due to detector status requirements. 

The number of events observed agrees with the number of events expected from e + e~ — > 7/177(7) 
plus the background. The two Monte Carlo generators give similar results although the KORALZ 
generator has a systematically lower cross-section than NNGG03/NUNUGPV. Following the discussion 
in section 2, the KORALZ Monte Carlo sample is used for all subsequent measurements and results 
concerning the e + e _ — > 1/177(7) process unless explicitly mentioned otherwise. 

Systematic errors on the cross-section measurement arising from uncertainties on the electromag- 
netic calorimeter energy scale and resolution, the description of the detector material and consequent 
conversion probabilities of photons in the central detector volume and coil, the integrated luminos- 
ity measurement, and the detector occupancy estimate, have been considered and evaluated to be 
negligible with respect to the statistical error. A relative systematic error of 4% is assigned to the 
cross-section measurement. This uncertainty comes dominantly from the estimated uncertainty on the 
efficiency based on comparing the different event generators. The modelling of the e + e~ — > 1/2/7(7) 
event fraction with a second photon detected in the forward detectors (| cos#| > 0.984), and therefore 
rejected by the forward energy vetoes is expected to be the main reason for the observed efficiency 
differences reported in Table 1. The cross-sections as a function of centre-of-mass energy are plotted 
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in Figure ^. The curves show the predicted cross-sections from the KORALZ event generator and 
the NUNUGPV analytical calculation for the Standard Model process e + e~ — ► 1/177(7). The data are 
generally consistent with the predictions but do not favour either estimate. 

In Figure ||a, the scaled energy of the most energetic photon is plotted against the cosine of its 
polar angle for events in the yfs=l72 GeV sample. The data are distributed as expected from the 
e + e~ — > ^1/7(7) Monte Carlo. Similar agreement is seen for the 130, 136 and 161 GeV data. In 
Figure ^b the polar angle distribution for the entire y / s=130-172 GeV sample is shown and agrees 
with the e + e~ — > 1/177(7) Monte Carlo expectation. If one calculates the mass recoiling against the 
photon (or against the two-photon system) in these events, one expects a peak in the recoil mass at 
Mz, since the vv predominantly comes from the decay of a Z°. One clearly sees this feature in the 
data as shown in Figure |j. In general there is also good agreement between data and Monte Carlo in 
this distribution. However, at \fs = 136 GeV, there is one event with a measured photon energy of 
84 GeV. Its estimated recoil mass is imaginary as the measured energy exceeds the beam energy, and 
it is shown in the distribution as occuring at zero recoil mass. A careful study of this event shows 
strong evidence, besides the measured energy, that it comes from a cosmic ray, well out-of-time with 
respect to the LEP beam crossing. In fact, it is sufficiently out-of-time so as to miss detection in 
several OPAL detector elements including the TOF system. It is left in the data sample because it 
passes all the selection criteria. It is not, however, considered to be a good physics event candidate. 

The single-photon selection was designed to allow for the presence of a second photon in order 
to accept events from the e + e~ — ► ^77(7) process. In the data 12 out of the 138 selected events 
are considered to be two photon events (i.e., have a second photon with deposited energy exceeding 
300 MeV). This is consistent with the expectation of 7.1 events from the KORALZ Monte Carlo. Ten 
of the 12 events are in common with the acoplanar-photons event selection. 

4.1.1 Search for e+e - — > XY, X — »■ Y7 - General case: M Y > 

Selected events at a given centre-of-mass energy are classified as consistent with a given value of Mx 
and My if the energy of the most energetic photon falls within the region kinematically accessible 
to a photon from the process e + e~ — * XY, X — ► Y7. The kinematic consistency criterion includes 
allowance for the energy resolution and incurs an inefficiency of less than 2% for all values of X and 
Y masses while accepting only those Standard Model e + e~ — > 2/177(7) events that are kinematically 
consistent with a given X and Y mass hypothesis. 

The kinematic region with true recoil mass significantly below Mz has only a small background 
expectation from Standard Model e + e~ — ► 1/177(7), but it may be populated as a result of energy mis- 
measurement in regions of the detector with poorer energy resolution. We reject events as candidates 
for XY production if the most energetic photon is in one of the following angular regions, 0.785 < 
I cos 9 I < 0.815 and | cos# | > 0.94, and the recoil mass is below 75 GeV. One data event is rejected 
by this cut compared with 1.4 events expected from the Standard Model e + e~ — > 1/1/7(7) process. We 
have studied the effect of further cuts to reduce the Standard Model e + e~ — > 1/177(7) background. 
We find that a significant improvement in the expected sensitivity can be achieved for small My by 
accepting only events with recoil mass significantly below Mz- Events are retained as candidates for 
small My, defined as My < 14+0. lMx (GeV), if the measured recoil mass is less than 75 GeV. For the 
complementary, large My region, My > 14 + 0.1Mx (GeV), no other cuts are applied. The boundary 
between the small and large My region was chosen so as to optimise the expected sensitivity^ for the 

3 The optimisation condition chosen was that the upper limit that one would expect to set, on average, in the absence 
of new physics contributions should be minimised. This definition has the advantage that it does not require one to 
specify the cross-section of possible new physics. 
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combined data sample. For simplicity, the same boundary is applied at each centre-of-mass energy. 

The selection efficiencies at each generated grid point for the SUSYGEN Monte Carlo events are 
shown for \fs = 130 GeV in Tableland for yfs = 172 GeV in Table [||. The efficiencies at intermediate 
centre-of-mass energies lie between those shown for these two centre-of-mass energies. These values 
include the efficiency of the kinematic consistency selection criteria which is higher than 98% at each 
generated mass point. The number of selected events consistent with each (Mx,My) value is shown 
in Figure ||] and can be compared with the number expected from Standard Model e + e~ — > 1/1/7(7) 
events as shown in Figure o. The background event described earlier does not survive the kinematic 
consistency cuts for any kinematically accessible point in the (Mx, My) plane. In general there is good 
agreement, and we proceed to set upper limits at 95% CL on the cross-section times branching ratio 
a(e + e~ — > XY)-BR(X — ► Y7). These upper limits are first calculated separately at each centre-of-mass 
energy, as shown in Figure |7[ A combined upper limit on the cross-section times branching ratio at yfs 
of 172 GeV is calculated combining the information from each centre-of-mass energy. The combination 
is performed assuming a cross-section centre-of-mass energy dependence of Px/s, where 0x is the speed 
of particle X. This combined limit is shown in Figure ||. The upper limits are calculated taking into 
account the expected number of Standard Model e + e~ — > 1/1/7(7) background events estimated from 
KORALZ using the method described in [|39]]. The estimated non-physics background is intentionally 
not taken into account in the limit calculation. The resulting combined upper limits range from 0.31 
pb to 1.8 pb. 

Systematic errors are due primarily to limited Monte Carlo statistics at the generated (Mx,My) 
points and the uncertainty on the efficiency parametrization across the (Mx, My) plane. The combined 
relative uncertainty on the efficiency is 4%. The effect of this uncertainty on the upper limits is 
calculated according to Q] and is found to be negligible. 



4.1.2 Search for e+e - — > XY, X — »■ Y7 - Special case: M Y « 

The above results include the case of My ~ that lies within the high mass-difference region in which 
the expected number of events is small. For example, for Mx = 100 GeV, two events are observed 
compared with an expected contribution from 1^7(7) of 0.8±0.1 events. For Mx = 170 GeV, no events 
are observed; the background expectation is 0.10 ±0.03 events. In this region, the requirement that the 
recoil-mass be less than 75 GeV eliminates all sources of physics background except a small residual 
contribution from e + e~ — ► 77(7); the expected contribution is about 0.01 events and is neglected. The 
upper limits for the My ~ function of Mx, are shown in Figure U. The resulting combined 

upper limits on cr(e + e - — > XY) • BR(X — > Y7) range from 0.36 pb to 0.76 pb. Interpretation of these 



results for the production of excited neutrinos is described in a separate paper [36|. 



4.2 Acoplanar-Photons 77(7) + $t 

After applying the acoplanar-photons selection criteria to the combined data sample, a total of 11 
events are observed. The predictions for the total number of Standard Model e + e~ — > 1/1/77(7) events 
is 6.3 ± 0.2, based on the KORALZ generator. Non-physics background as well as contributions from 
other Standard Model processes are negligible. The breakdown of the observed and expected number 
of events for the different centre-of-mass energies is given in Table |4|. Within the kinematic acceptance, 
the selection efficiency^ for e + e~ — > 1/1/77(7) events is rather constant as a function of centre-of-mass 

4 Before accounting for efficiency losses due to detector occupancy and detector status requirements. 
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energy, varying from about (67-70)% for the KORALZ generator and from (67-72)% for NNGG03 and 
NUNUGPV. The mean efficiencies are (68.6 ± 1.5)% for KORALZ and (70.4 ±0.5)% for NNGG03 and 
NUNUGPV. As the cross-section measurements are statistics limited, an efficiency of (69 ±3)% is used 
independent of y/s. Additional systematic errors arise due to the energy scale for low-energy photons 
(5%), and from uncertainty on the luminosity measurement (< 1%). The measured cross-sections for 
the process e + e~ — > 1/7777(7), within the kinematic acceptance defined by the energy and polar angle 
selection criteria described earlier, are included in Table |] as are the cross-section predictions from 
the KORALZ generator. 



There were no selected events having iV 7 > 2, compared with an expectation from KORALZ of 
0.34 ±0.06 events. The kinematic properties of the selected events are summarized in Table |5| and dis- 
played in Figures 0-13 where they are compared with the predicted distributions for e + e~ — > ^7/77(7), 



obtained using the KORALZ generator, normalized to the corresponding integrated luminosity. Fig- 
ure shows the recoil-mass distributions of the selected acoplanar-photon pairs. These are peaked 
near the mass of the Z° as expected for contributions from e + e~ — > 2/7777(7). The resolution of the 
recoil mass is typically about 2-3 GeV at each centre-of-mass energy for M recoi i « Mz- Figure p] 
shows the X2 vs. x\ distributions for each centre-of-mass energy and for the combined data sample. 
The projections of the scaled energy of the least-energetic photon are given in in Figure 12. Figure ^ 



shows the distributions of the invariant mass, M 77 , for the selected acoplanar-photon pairs at each 
centre-of-mass energy. The mass resolution is typically 0.6-1.4 GeV. A search for H° — ► 77 has been 



recently published by OPAL [17|. 



For the data at y/s = 161, 172 GeV, the measured distributions agree with the 2/1/77(7) expectation. 
For the y/s = 130, 136 GeV data there is an apparent excess of events. However, with the exception of 



the X2 distribution in Figure 12, for which the measured distribution peaks more strongly than expected 
at low values, the distributions appear consistent in shape with the expectation from e + e~ — > 1/7/77(7). 
This excess was also remarked on in a previous publication pj. 



4.2.1 Search for e+e - — > XX, X — »• Y7 - General case: M Y > 



Selected events are classified as consistent with a given value of Mx and My if each of the two selected 
photons falls within the region kinematically accessible to photons from the process e + e~ — ► XX, 
X — * Y7. As before, this includes allowance for resolution effects. Monte Carlo events were generated 
at each centre-of-mass energy. The selection efficiencies at each generated grid point for the y/s = 172 
GeV SUSYGEN Monte Carlo events are shown in Table ||. These values include the efficiency of 
the kinematic consistency selection criteria which is higher than 95% at each generated point in the 
(Mx,My) plane. Similar efficiencies are obtained at the other centre-of-mass energies. 

Figures || (a)-(d) show the 95% CL cr(e + e~ -> XX) • BR 2 (X -> Y7) exclusion plots obtained at 
y/s = 130, 136, 161 and 172 GeV respectively. Because of the current uncertainties on the modelling 
of the Standard Model background, as discussed earlier, these limits and the limits presented below for 
this topology, have been calculated without taking into account the background estimate. Events from 
e + e~ — > 1/7/77(7) are typically characterized by a high-energy photon from the radiative return to the 
Z° and a second lower energy photon. The kinematic consistency requirements, however, require that 
the two photons have energies within the same (kinematically accessible) region. For this reason, two 
of the selected events are inconsistent with any (Mx, My) point for Mx > 45 GeV and Mx — My > 5 
GeV. Figure 15 shows the 95% CL cr(e + e~ — ► XX) • BR 2 (X — * Y7) exclusion plot obtained from the 



combined data sample assuming that <r(e + e — > XX) scales with centre-of-mass energy as Px/s. For 
the combined plot, the maximum value of the limit is 0.80 pb. The minimum value is 0.18 pb. 
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Systematic errors are due primarily to limited Monte Carlo statistics at the generated (Mx, My) 
points and the uncertainty on the efficiency parametrization across the (Mx,My) plane. The com- 
bined relative uncertainty on the efficiency varies from about (3-6)% across the plane. All systematic 
uncertainties are accounted for in the manner advocated in reference [ p£p . This also applies to the 
limits for the My ~ case, presented in the next section. 



4.2.2 Search for e+e - — >• XX, X — > Y7 - Special case: My ~ 



For the special case of My ~ the kinematic consistency cuts applied differ from those used for the 



general case. One can calculate [14] the maximum and minimum masses, M^ ax and M™ n , which are 
consistent with the kinematic properties of the two photons, assuming a massless Y. As this argument 
is based only on kinematics it applies generally to the case where the acoplanar photon pair originates 
from pair production of heavy neutral particles which decay radiatively to massless invisible final 
states; e + e~ — > XX followed by X — > Y7, My ~ 0. These maximum and minimum mass values can 
provide further suppression of the ^77(7) background while retaining high efficiency for the signal 
hypothesis. The background suppression achieved with kinematic consistency requirements based on 
this procedure is much better than that obtained for the general case since, in this case, the full 
event kinematics are used. Figure [l6] shows M™ n vs. M-^ 1 ax for events passing the selection criteria 
described in section |3]^ for, (a) 1/1777(7) Monte Carlo and (b-d) e + e~ — ► XX, X — > Y7, My 
signal Monte Carlo, at -y/i = 172 GeV, for three values of Mx- The signal Monte Carlo distributions 
are dominantly populated at maximum mass values greater than or equal to the generated mass of X 
(e.g. v* or Xi)- Resolution effects shift some entries to lower masses. Requiring that the maximum 
kinematically allowed mass be greater than Mx — 5 GeV retains more than 95% relative efficiency for 
signal at all values of Mx while suppressing 41% to 96% of the remaining ^1777(7) events. Similar 
efficiencies are obtained at each of the other centre-of-mass energies. 

The kinematic properties of the selected events, shown in Table |5], include the values of My^**. The 



data distributions in M™ 11 vs. M™ are shown for each of the centre-of-mass energies in Figure 17. In 
each case the distribution for ^7777(7) events is also shown. For illustrative purposes, the efficiencies 
calculated from Monte Carlo events generated at 172 GeV are shown in Table before and after 
application of the cut on M^ iax . Also shown is the ^7777(7) rejection efficiency obtained with the 
Mx ax cut. Signal efficiencies at other centre-of-mass energies are similar. 

Based on the efficiencies and the number of events observed at each centre-of-mass energy, we 
calculate 95% CL upper limits on c(e +e _ -» XX) • BR 2 (X -> Y 7 ) (for My « 0) clS £1 function of Mx- 



These are shown in Figure 18. Also shown is the 95% CL upper limit obtained from the combined data 
sample, assuming a centre-of-mass energy dependence of the cross-section of /3x/s. This combined 
limit is 0.5 pb or less for values of Mx from 45 GeV up to the kinematic limit. 

To set combined, model dependent limits on the mass of the Xi (NLSP) in supersymmetric models 
in which the LSP is a light gravitino, we sum the number of observed events consistent with each value 



of Mx- This distribution is shown in Figure IS where the solid line shows the number of observed 
candidates consistent with a given value of Mx and the dashed-dotted line shows the background 
expectation from KORALZ. The number of candidate events is consistent with the number of back- 
ground events expected from 1/1777(7). The thick solid line shows the 95% confidence level upper limit 
at each mass value. The background expectation is not taken into account when calculating the limit. 
Also shown in Figure [HI are the numbers of events expected from the Lopez and Nanopoulos no-scale 
supergravity model [O] and from the model of Babu, Kolda and Wilczek in which the neutralino 
composition is purely gaugino (bino). For both of these models, the cross-section has been evaluated 
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at Born level. Based on these distributions, these two models are ruled out at the 95% confidence 
level for M^o < 61.3 GeV and 69.4 GeV, respectively. 

As described above, the efficiencies over the full angular range have been calculated using isotropic 
angular distributions for production and decay of X. The validity of this model has been examined 



based on the angular distributions calculated for photino pair production in [12|. For models proposed 
in [ 13 1 , the production angular distributions are more central and so this procedure is conservative. 
For a 1 + cos 2 9 production angular distribution, expected for t-channel exchange of a very heavy 
particle according to [p^| , the relative efficiency reduction would be less than 2% for all points in the 
M x , M Y plane (for M x - M Y > 5 GeV). 

Interpretation of these results for the production of excited neutrinos is described in a separate 
paper [36]. 



5 Conclusions 



We have searched for photonic events with large missing energy in two different and complementary 
topologies in data taken with the OPAL detector at LEP, at centre-of-mass energies in the region of 
130-172 GeV. 

In the single-photon selection which requires at least one photon with xt > 0.05 in the region 
| cos 9 | < 0.82 or one photon with xt > 0.1 in the region 15° < 9 < 165° (|cos#| < 0.966) a total 
of 138 events are observed in the data compared to the KORALZ prediction for the contribution 
from the Standard Model process e + e~ — * vV r y('y) of 141.1 ± 1.1 events and an expected non-physics 
background of 2.3 ±1.1 events. The corresponding cross-sections for e + e~ — ► ^7(7) are 10.0 ± 2.3, 
16.3 ± 2.8, 5.3 ± 0.8 and 5.5 ± 0.8 pb for yfs = 130, 136, 161 and 172 GeV, respectively, in agreement 
with the Standard Model expectations. We derive upper limits on the cross-section times branching 
ratio for the process e + e~ — > XY, X —* Y7 for the general case of massive X and Y. The limits vary 
from 0.31 to 1.8 pb in the region of interest of the (Mx,My) plane and include the special case of 
My ~ 0, where the limit varies between 0.36 and 0.76 pb for the M x mass range from Mz to 172 
GeV. 

The acoplanar-photons selection requires at least two photons with scaled energy x 7 > 0.05 within 
the polar angle region 15° < 9 < 165° or at least two photons with energy E~ > 1.75 GeV with one 
satisfying | cos 9 \ < 0.8 and the other satisfying 15° < 9 < 165°. From the combined data sample 11 
events are selected. The KORALZ prediction for the number of events from e + e~ — ► 1^77(7) is 6.3 ± 
0.2. The cross-section for this process is measured at each centre-of-mass energy (see Table |). Due 
to the uncertainties in the current modelling of the Standard Model background, e + e~ — > 1/1/77(7), an 
limits from the acoplanar-photons analysis were calculated without taking into account the background 
estimate. Based on a variety of kinematic distributions, all observed events appear consistent with 
e + e~ — > 1/1777(7). We derive 95% CL upper limits on a (e+e - -» XX) • BR 2 (X -> Y7) ranging from 
0.18 to 0.80 pb for the general case of massive X and Y, and from 0.35 to 0.50 pb for the special case 
of M Y 0. 

For the single-photon and acoplanar-photons search topologies, the general case of massive X and 
Y is relevant to the supersymmetry models in which X = X2 an d V = Xii with x\ ~^ Xil an d X? 
stable. The special case of My ~ is of particular interest for single and pair production of excited 
neutrinos and for supersymmetric models in which the LSP is a light gravitino and the Xi is the NLSP 
which decays as x? — * G7. For the latter scenario, the results of the acoplanar-photons search are 
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used to place model-dependent lower limits on the Xi mass. Comparison with the model predictions 
of Lopez and Nanopoulos [14| permits exclusion of that model for M^o < 61.3 GeV. A similar model 
from Babu, Kolda and Wilczek [16| is excluded for M^o < 69.4 GeV. The results of these searches 
have also been used to place limits on the production of excited neutrinos 
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x 11 DKg 


e K \ /o > 
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19 


25.3 ± 0.4 


27.5 ± 0.3 


0.2 ± 0.2 


81.6 ± 0.6 


79.6 ± 0.5 


10.0 ± 2.3 


136 


34 


23.3 ± 0.4 


25.8 ± 0.3 


0.3 ± 0.2 


79.7 ± 0.7 


80.0 ± 0.5 


16.3 ± 2.8 


161 


40 


48.3 ± 0.6 


51.3 ± 0.4 


0.9 ± 0.8 


75.2 ± 0.5 


72.3 ± 0.4 


5.3 ± 0.8 


172 


45 


44.3 ± 0.6 


46.1 ± 0.4 


0.9 ± 0.8 


77.9 ± 0.5 


74.5 ± 0.4 


5.5 ± 0.8 


130-172 


138 


141.2 ± 1.1 


150.7 ± 0.7 


2.3 ± 1.1 





Table 1: For each centre-of-mass energy, the table shows the number of events observed in the OPAL 
data, the number expected based on the KORALZ (K) and NUNUGPV/NNGG03 (TV) e+e~ -► 1/1/7(7) 
event generators and the number of events expected from non-physics backgrounds. Also shown are the 
efficiencies obtained from the two generators, within the kinematic acceptance of the single-photon 
selection, and the measured cross-sections within the kinematic acceptance, determined using the 
efficiencies obtained with the KORALZ generator. The quoted cross-section errors are statistical. 



M x + M Y 
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M Y = M x /2 


M Y = M x - 15 GeV 


M Y = M x - 5 GeV 


110 


80.2 ± 1.7 


83.7 ± 1.5 


80.9 ± 1.7 


51.3 ± 2.2 


120 


82.3 ± 1.6 


86.6 ± 1.4 


80.2 ± 1.7 


49.2 ± 2.2 


128 


84.9 ± 1.5 


84.4 ± 1.5 


82.5 ± 1.6 


49.6 ± 2.2 



Table 2: Single photon selection efficiency (%) as a function of the sum of Mx and My versus various 
My values for the process e + e~ — > XY, X — > Y7. These efficiencies are for yfs = 130 GeV. Masses 
given are in GeV. 
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M Y = M x - 5 GeV 
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75.3 ± 1.3 


32.2 ± 1.5 
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86.5 ± 1.0 


74.0 ± 1.3 


30.9 ± 1.4 


170 


82.9 ± 1.1 


88.7 ± 0.9 


74.8 ± 1.3 


30.8 ± 1.4 



Table 3: Single photon selection efficiency (%) as a function of the sum of Mx and My versus various 
My values for the process e + e~ — > XY, X — > Y7. These efficiencies are for y/s = 172 GeV. Masses 
given are in GeV. The errors are statistical. 
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0.33 ± 0.02 
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6.29 ± 0.24 





Table 4: Number of acoplanar-photons events observed and expected at each centre-of-mass energy 
region and the corresponding measured and expected (KORALZ) cross-sections within the kinematic 
acceptance of the acoplanar-photons event selection. 
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Table 5: Kinematic properties of the events passing the acoplanar-photons selection. All masses are 
given m GeV. Mf ax is defined in Section gjjg . 
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Table 6: Acoplanar-photons selection efficiencies (%) for the process e + e~ — > XX, X — > Y7 at y/s = 172 
GeV for various Mx and My (in GeV) . These values include the efficiency of the kinematic consistency 
cuts. Similar efficiencies are obtained at each value of y/s. The errors are statistical. 
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Table 7: Acoplanar-photons event selection efficiency (%), as a function of mass, for the process 
e + e~ — * XX, X — ► Y7, for My ~ 0. These numbers are for y/s = 172 GeV. The first column shows the 
efficiency of the selection described in section 3.2. The second column shows the selection efficiency 
after the cut on M™ described in section 4.2.2. The final column shows the rejection efficiency (%) 
of the Mjp 8 * > Mx — 5 GeV cut for ^77(7) events. The errors are statistical. 
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Figure 1: For single-photon conversion candidates, a) shows the maximum angular separation in 9 
and 4> of the primary photon candidate and the highest px track in the event. For the non-conversion 
single-photon selection, b) shows the difference between the observed TOF timing and the timing 
expected for a photon from the interaction point, for data events passing all cuts or failing only the 
timing cut. For the acoplanar-photons selection, c) shows the maximum cluster extent for data events 
failing only the anti-77(7) cuts (shaded histogram) and for data events failing the TOF requirements 
or the TOF requirement and the special background vetoes (unshaded histogram). For the acoplanar- 
photons selection, d) shows the distribution of the acoplanarity angle for data events failing only the 
total energy cut and/or the cut on pt{h) (predominantly e + e~ — > 77). In a) and d) the solid points 
with error bars show the OPAL data while the overlaid histograms represent the expectation from 
2477(7) and e + e~ — > 77(7) Monte Carlo, respectively, normalized to the luminosity of the data. 



22 



1 — I — i — i — i — i — I — i — i — i — r 



1 — i — i — i — I — i — i — i — i — I — i — r 



OPAL 



1 — i — i — r 



e + e — > y (y)+invisible particle(s) 



NUNUGPV prediction 
KORALZ prediction 



J I I L 



J I I L 



J I I L 



J I I L 



J L 



1 



J I I L 



120 130 140 150 



160 170 180 190 200 
Centre-of-mass energy (GeV) 



Figure 2: The measured value of <r(e + e - — > 7(7) + invisible particle(s)), within the kinematic 
acceptance of the single-photon selection, as a function of yfs. The data points with error bars are 
OPAL measurements at centre-of-mass energies of 130, 136, 161 and 172 GeV. The curves are the 
predictions for the Standard Model process e + e~ — > 1^7(7) from the KORALZ generator and the 
NUNUGPV analytical calculation. 
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Figure 3: a) Distribution of x 7 vs cos 6 for the most energetic photon in the single photon selection. 
The fine points are the KORALZ e + e~ — > 2/177(7) Monte Carlo and the solid triangles are the data. 
This plot is for y/s = 172 GeV. b) The cos 9 distribution for the most energetic photon in the single 
photon selection. The points with error bars are the data and the histogram is the expectation from 
the KORALZ e + e~ — > vV^(^) Monte Carlo normalized to the integrated luminosity of the data. This 
plot is for the combined data set \Js = 130-172 GeV. 



24 



OPAL 



> 
O 

1/1 c 

3 5 







~i — i — i — | — i — i — i — | — i — i — i — | — i — i — i — | — i — i — i — | — i — i — i — | — i — i — i — | — i — i — i — | — i — i — i — | — i — i — r 

Vs=130-136 GeV 



L 



J I I I I I I I I I I L 




1 I I 



j i i i_ 



20 40 60 80 100 120 140 160 180 200 

Recoil Mass (GeV) 



% 8 

"J 6 
<s 

^ 4 

+^ 

a 

> 

H 



~ i — i — i — i — i — i — i — i — i — i — i — i — i — i — i — i — r 

Vs=161GeV 



_l I I I I I I I I I I I L 




i — i — | — i — i — i — | — i — i — i — | — i — i — i — | — i — i — i — | — i — i — r 



T 



T 



_l I I I I I l_ 



20 40 60 80 100 120 140 160 180 200 

Recoil Mass (GeV) 



> 

a 6 

» 4 

a ~ 

o> 2 

W 



^30 
O 

">20 

^ 10 

> 

H 



~i — i — i — i — i — i — i — i — i — i — i — i — i — i — i — i — i — i — i — i — i — i — i — i — i — i — i — i — i — i — i — i — i — i — i — i — r 

Vs=172 GeV 



J I I I I I I I I L=J 




20 40 60 80 100 120 140 160 180 200 

Recoil Mass (GeV) 



—1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 

Vs=130-172 GeV 




1 1 1 — 


; i J 

\ i i i I i i i I i i i i i_<4-'-r» 


V < 


1 1 1 



20 40 60 80 100 120 140 160 180 200 

Recoil Mass (GeV) 



Figure 4: The recoil mass distribution for events passing the single photon selection for the y/s = 
130-136 GeV, 161 GeV, 172 GeV, and combined 130 - 172 GeV data samples. The points with error 
bars are the data and the histograms are the expectations from the KORALZ e + e~ — > 1/1/7(7) Monte 
Carlo normalized to the integrated luminosity of the data. 
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Figure 5: Number of single photon candidate events in the combined data sample (^/s = 130-172 
GeV) consistent with each set of mass values (M x , My) for the process e + e~ — > XY, X — > Y7, after 
application of all selection criteria including kinematic consistency requirements. Lines are drawn 
around the boundaries defined by M x + My = 172 GeV, M x = My, and M x + My = M z , and to 
display the boundary between the small and large My regions. Regions defined by M x — My < 5 
GeV and M x + My < Mz are not considered in limit calculations, and no events are displayed in 
these regions. 
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M x (GeV) 

Figure 6: Number of Standard Model ^7(7) events predicted by KORALZ to pass all single photon 
selection criteria for the combined data sample (y/s = 130-172 GeV), which are consistent with the 
process e + e~ — > XY, X — ► Y7 at each set of mass values (Mx, My). This figure gives the expected 
Standard Model contribution to Figure [|. Boundaries and delineated regions are the same as in 
Figure [|. 
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Figure 7: 95% CL upper limits on <r(e + e _ — ► XY) • BR(X — > Y7) as a function of Mx and My at 
each centre-of-mass energy. Lines are drawn to indicate the kinematically allowed boundaries defined 
by Mx + My = yfs. Other boundaries and delineated regions are as defined for Figure [|. 
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Figure 8: The 95% CL upper limit on cr(e+e~ -> XY) • BR(X -> Y7) at ^/s = 172 GeV as a function 
of Mx and My, obtained from the combined data sample assuming a cross-section scaling of fix/s. 
The boundaries and delineated regions are as defined for Figure || 
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M x (GeV) 

Figure 9: The 95% CL upper limits on <r(e + e - — ► XY) • BR(X — ► Y7) as a function of Mx, assuming 
My ~ 0, for each centre-of-mass energy. The upper limit on a(e + e~ — > XY) • BR(X — > Y7), evaluated 
at -y/s = 172 GeV, obtained from the combined data sample is also shown. The combination was 
performed assuming a cross-section scaling of (3x/s. The allowed region is shaded. 
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Figure 10: Recoil-mass distributions for the selected acoplanar-photons events for each centre-of-mass 
region and for the combined data sample. The data points with error bars represent the selected OPAL 
data events. In each case the shaded histogram shows the expected contribution from e + e~ — > 2/1777(7) 
events, normalized to the total integrated luminosity. The KORALZ generator was used for these 
distributions. 
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Figure 11: X2 vs x\ for the selected acoplanar-photons events for each centre-of-mass region and for the 
combined data sample. The large points represent the selected OPAL data events. The smaller dots 
show the expected distribution for events from the process e + e~ — > 1/1/77(7) events, from KORALZ. 
The normalization of the Monte Carlo distributions is arbitrary. However, for the combined plot the 
relative normalizations of the various data samples is maintained. 
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Figure 12: Distributions of scaled energy of the lower energy photon for the selected acoplanar- 
photons events for each centre-of-mass region and for the combined data sample. The data points 
with error bars represent the selected OPAL data events. In each case the shaded histogram shows the 
expected contribution from e + e~ — > ^77(7) events, normalized to the total integrated luminosity. 
The KORALZ generator was used for these distributions. 
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Figure 13: Distributions of the 77 invariant mass of the accepted acoplanar-photon pairs for each 
centre-of-mass region. The points with error bars represent the OPAL data while the shaded his- 
tograms show the predicted distributions for events from e + e~ — > 1/1/77(7) events, from KORALZ, 
normalized to the corresponding integrated luminosity. 
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Figure 14: For each centre-of-mass energy, the shaded areas show 95% CL exclusion regions for 
cr(e+e- -► XX) • BR 2 (X -► Y 7 ). No limit is set for mass-difference values Mx — My < 5 GeV, defined 
by the lower line above the shaded regions. The upper line is for Mx = My. 
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Figure 15: The shaded areas show 95% CL exclusion regions for <r(e+e- -► XX) • BR 2 (X -► Y 7 ) at 
\fs = 172 GeV, obtained from the combined data sample assuming a cross-section scaling of flx/s. 
No limit is set for mass-difference values Mx — My < 5, defined by the lower line above the shaded 
regions. The upper line is for Mx = My. 
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Figure 16: M£ lin vs. Mf ax for (a) 1/1777(7) Monte Carlo (KORALZ) and (b)-(d) e+e" -► XX, 
X — > Y7, My ~ Monte Carlo (SUSYGEN) for various Mx- The vertical lines represent the chosen 
cut values discussed in the text. These are for Monte Carlo generated at a centre-of-mass energy of 
172 GeV. 
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Figure 17: M™ 11 vs. M™ ax for the accepted acoplanar-photons events for each centre-of-mass energy 
region. Overlaid are the expected distributions for contributions from the Standard Model process 
e + e~ — ► ^1/77(7), from the KORALZ generator. 
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Figure 18: The 95% confidence level upper limit on <r(e + e - — ► XX) • BR 2 (X — > Y7) for the case 
My ~ 0, as a function of Mx, for each value of sjs. Also shown is the limit, evaluated at y/s = 172 
GeV, obtained from the combined data sample assuming a cross-section scaling of (3x/s. 
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Figure 19: Results for e + e~ — > XX, X —* Y7, with My ~ for the combined data sample. The lower 
solid line shows the number of observed candidates consistent with a given value of Mx- The dashed- 
dotted line shows the expected contribution from 1/7777(7) obtained from the KORALZ generator. The 
thick solid line shows the 95% confidence level upper limit on the number of candidate events. The 
dashed (dotted) line shows the expected number of events from the model of Lopez and Nanopoulos []l4|] 
(Babu, Kolda and Wilczek |jl§). Within these models, Xi masses less than 61.3 (69.4) GeV are 
excluded at the 95% confidence level. 
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